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Abstract The unitary representation of exact quantum Poincare mapping is constructed. 
It is equivalent to the compact representation | jTT| [22] , in a sense that it yields equiva- 
lent quantization condition with important advantage over the compact version: since it 
preserves the probability it can be literally interpreted as the quantum Poincare mapping 
which generates quantum time evolution at fixed energy between two successive crossings 
with surface of section (SOS). SOS coherent state representation (SOS Husimi distribu- 
tion) of arbitrary (either stationary or evolving) quantum SOS state (vector from the 
Hilbert space over the configurational SOS) is introduced. Dynamical properties of SOS 
states can be quantitatively studied in terms of the so called localization areas which 
are defined through information entropies of their SOS coherent state representations. 
In the second part of the paper I report on results of extensive numerical application 
of quantum SOS method in a generic but simple 2-dim Hamiltonian system, namely 
semiseparable oscillator. I have calculated the stretch of 13 500 consecutive eigenstates 
with the largest sequential quantum number around 18 million and obtained the follow- 
ing results: (i) the validity of the semiclassical Berry-Robnik formula for level spacing 
statistics was confirmed and using the concept of localization area the states were quanti- 
tatively classified as regular or chaotic, (ii) the classical and quantum Poincare evolution 
were performed and compared, and expected agreement was found, (iii) I studied few 
examples of wavefunctions and particularly, SOS coherent state representation of regu- 
lar and chaotic eigenstates and analyzed statistical properties of their zeros which were 
shown on the chaotic component of 2-dim SOS to be uniformly distributed with the 
cubic repulsion between nearest neighbours. 
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1 Introduction 



In recent papers |L9|, |22|, the author established the surface of section (SOS) reduc- 
tion of quantum dynamics in complete analogy with the famous Poincare SOS reduction 
of classical dynamics ||16|| . In a similar way, as the smooth bound autonomous classical 



Hamiltonian dynamics in 2/-dim phase space can be reduced to a discrete area preserving 
(Poincare) mapping over 2(f — l)-dim SOS, the essential information about quantum dy- 
namics is contained in the energy-dependent propagator which acts over the Hilbert space 
of complex valued functions over (/ — l)-dim configurational surface of section (CSOS) 



- the so called quantum Poincare mapping (QPM). One can prove [|T^, |22|, |23|] the fol- 
lowing interesting results which constitute quantum SOS reduction: (i) the eigenenergies 
of the original bound Hamiltonian are the points where the QPM possesses "nontrivial" 
fixed points, (ii) there exist additional three energy-dependent propagators which prop- 
agate from/to L 2 (CS) to/from L 2 (CSOS) or L 2 (CS) (interpreted as propagators at fixed 
energy and without crossing the CSOS) and decompose the energy-dependent quantum 
propagator (E — H)^ 1 in terms of newly defined propagators with respect to CSOS. 
One can consistently write the QPM as the product of two generalized (non-unitary) 
multi-channel scattering operators which correspond to two scattering problems which 
are obtained by cutting one half of CS along CSOS off and attaching a semi-infinite flat 
waveguide instead, in such a way that the two scattering Hamiltonians remain contin- 
uous. Such compact QPM (since, in a language of functional analysis, it is a compact 
operator) is consistent with its semiclassical limit |23f derived by Bogomolny a few years 
ago ||. The same form of QPM was proposed for exact quantization of billiards by 
Smilansky and coworkers ]3T], . And quite recently, Rouvinez and Smilansky []30| used 
the same technique for the quantization of a smooth 2-dim Hamiltonian. 

But this compact QPM has one important disadvantage: Since it is not norm- 
preserving i.e. nonunitary operator it cannot be used to define quantum SOS evolution 
in analogy with the classical time evolution generated by the classical Poincare map- 
ping. It is only approximately unitary (its eigenvalues typically lie in the neighbourhood 
of a complex unit circle) and becomes unitary only in the semiclassical limit ||. It 
is the purpose of this paper to show that there exists another consistent, equivalent 
and unique quantum SOS reduction in which QPM is strictly unitary. Such QPM is 
norm-preserving and it may be used to define unique quantum SOS time evolution of 
quantum SOS states (vectors from L 2 (CSOS)) which corresponds to classical SOS time 
evolution. Each quantum SOS state may be represented by a kind of quantum SOS 
probability distribution function which is defined as its coherent state representation - 
Husimi distribution in terms of suitably defined SOS coherent states. This SOS Husimi 
distribution (analogously one could define nonpositive SOS Wigner distribution) may be 
used to study localization properties of quantum SOS states and their correlation with 
classical invariant components of SOS. I define a quantitative measure of localization, 
the so called localization area of a quantum SOS state which is defined through the 
information entropy of the corresponding SOS Husimi distribution. 

The second purpose of this paper is to demonstrate the practical power of the uni- 
tary quantum SOS method as applied to a geometrically special, simple but dynamically 
generic (nonlinear) autonomous Hamiltonian system, namely 2-dim semiseparable oscil- 
lator. Using SOS reduction of quantum dynamics we were able to calculate a stretch 
of 13 500 consecutive eigenenergies and eigenstates with the sequential quantum number 
around 18 million. We have chosen a generic regime where classical dynamics of a system 
is mixed with regular and chaotic regions coexisting in phase space and on SOS. Thus I 
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was able to review many phenomena and confirm several conjectures of quantum chaos: 
(i) I confirmed the validity of semiclassical Berry- Robnik formula @] for the energy level 
spacing distribution and quantitatively classified the eigenstates as regular or chaotic 
by means of their localization areas, (ii) I have explicitly studied SOS time evolution 
of quantum and classical SOS probability distributions and found expected agreement 
within the so called break time, (iii) I have studied few special examples of wavefunctions 
of eigenstates in configuration space, (iv) and SOS Husimi distributions of eigenstates 
and analyzed the statistical properties of their zeros and found that they are uniformly 
distributed over the chaotic component of SOS with the cubic repulsion between nearest 
neighbours while in regular region they typically condense close to 1-dim classically in- 
variant or anti-Stokes [17| curves. This confirms and extends the conjecture by Leboeuf 
and Voros 



2 Unitary Quantum Surface of Section Method 



First we shall review some of the already known basic results [23] of the quantum SOS 
method which will be needed for further derivation of its unitarized version. We study 
autonomous and bound (at least in the energy region of our concern) Hamiltonian systems 
with few, say /, freedoms, living in an /—dim configuration space (CS) C. One should 
also provide a smooth orientable (/ — l)-dim submanifold of CS C which shall be called 
configurational surface of section (CSOS) and denoted by So. CSOS So cuts the CS 
C in two pieces which will be referred to as upper and lower and denoted by C a with 
two values of the binary index a =|, j. The two parts of the boundary of C which lie 
above/below So will be denoted by B a . If CS C is infinite then q 6 B a will stand for the 
limiting process \q\ — ► oo with q e C a . In arithmetic expressions the arrows will have 
the following values \— +1, {— —1. We choose the coordinates in CS, q = (x,y) G C in 
such a way that the CSOS is given by a simple constraint y = 0, or S = (S, 0). These 
coordinates need not be global, i.e. they need not uniquely cover the whole CS, but they 
should cover the open set which includes the whole CSOS So- This means that every 
point in So should be uniquely represented by CSOS coordinates x e S which may be 
more general than Euclidean coordinates TZ^~ l (e.g. (/ — 1) dim sphere S^ 1 ). 
We shall assume that the Hamiltonian has the following quite general form 

H = ^-p 2 y + H\p x ,x,y), (1) 

so that the kinetic energy is quadratic at least perpendicularly to the CSOS. There is 
straightforward generalization of the theory to the cases of nonconstant but positionally 
dependent mass and slightly less straightforward generalization to the cases of external 



gauge fields where we have also terms which are linear in momentum p y p3 . 

In quantum mechanics, the observables are represented by self-adjoint operators in a 
Hilbert space 7i of complex- valued functions ^{q) over the CS C which obey boundary 
conditions ^f(dC) = and have finite L 2 norm f c dq\^(q)\ 2 < oo. We shall use the Dirac's 
notation. Pure state of a physical system is represented by a vector — ket which 
can be expanded in a convenient complete set of basis vectors, e.g. position eigenvectors 
\q) = \x,y), = f c dq\q)(q\ty) = f c dq^/(q)\q) (in a symbolic sense, since \q) are 
not proper vectors, but such expansions are still meaningful iff ^(q) = (q\^) is square 

2 They exist since SOS Husimi function may be written according to Bargmann Q in terms of analytic 
functions of (/ — 1) complex variables. 
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integrable i.e. L 2 (C)-function). Every ket \^) G 7i has a corresponding vector from the 
dual Hilbert space 7i', that is bra (\&| e TC', (^\q) = (q\^f)*. We shall use mathematical 
accent "to denote linear operators over the Hilbert space 7i. An operator A acts either 
on abstract ket or on the corresponding wavefunction A^(q) = (q\A\ty). The major 
problem of bound quantum dynamics governed by the self-adjoint Hamiltonian H is to 
determine the eigenenergies E for which the Schrddinger equation 

H^ a (q,E) = E^ a (q,E) (2) 

has nontrivial normalizable solutions — eigenf unctions ty(q,E). 
Operators of SOS coordinates x and p x , defined by 

(x,y\x\^f) = xW(x,y), 
(x,y\p x \V) = -ihd x ^(x,y), 

can be viewed also as acting on functions ip{x) of x only and therefore operating in some 
other, much smaller Hilbert space of square-integrable complex-valued functions over a 
CSOSSo 

{cclac^} = XTp(x), 
{x\p x \ip} = -ihd x ip(x). 

Vectors in such reduced SOS Hilbert space, denoted by £, will be written as and linear 
operators over £ will wear mathematical accent "like restricted position x and momentum 
p x . Eigenvectors \x} of SOS position operator x provide a useful complete set of basis 
vectors of L. The quantum Hamiltonian can be written in position representation at 
least locally as 

H = -^-d 2 y +H / (y), H , (y) = H\- i hd x ,x,y). (3) 

The eigenstates of the reduced Hamiltonian H' '(0) = H'(0)\c restricted to the SOS Hilbert 
space £, \n) G £ 

F'(0)|n} = £>}, (4) 

which are called SOS eigenmodes, provide a useful (countable n = 1, 2, . . .) complete and 
orthogonal basis for C since H'(0) is a self-adjoint operator with discrete spectrum when 
its domain is restricted to C. 



2.1 Scattering formulation of quantum SOS 

In order to make the paper selfcontained we shall here review some crucial ingredients 



of a scattering formulation of quantum SOS method 0? H, £(J which will be needed 
for further derivation of its unitarized version. 

Let us consider for the moment an open quantum (quasi 1-dim) waveguide with the 
Hamiltonian 



H iree = i + H'(p x , x, 0) = -^-d 2 y + H'(0). (5) 
2m 2m y 



The motion it describes is free in y— direction and bound in all other ^-directions. The 
basic solutions of the Schrddinger equation at some arbitrary energy E in such waveg- 
uide are separable to products of plane waves in y— direction and SOS eigenmodes 
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in cc— direction {x\n}e ±lkn ^ v , with the corresponding wavenumber determined by the 
energy difference E — E' n available for the free motion perpendicular to the CSOS 

k n (E) = 

Such basic solutions for any n will be called channels, and there are typically finitely many 
open or propagating channels n for which the energy difference E — E' n is positive, that is 
the wavenumber k n (E) is positive real, and infinitely many closed or evanescent channels 
n for which the energy difference E—E' n is negative, which means that wavenumber k n (E) 
is positive imaginary. Sometimes I will write C (E) for the finitely dimensional subspace 
of C spanned by SOS eigenmodes corresponding to open channels and C C (E) for its 
orthogonal complement, so that £ = C Q (E) © C C (E). If one defines the wavenumber 
operator over C 




K(E) = £ k n (E)\n}{n\ = J -^(E - H'(0)). (6) 



then a general (nonbound) solution of the Schrodinger equation in the waveguide may 
be compactly written in terms of arbitrary two SOS states G C 

V(x, y, E) = - 7 ^={x\K- 1/2 (E) \e lR{E)y \^} + e^*^^}] (7) 



■im 



where the prefactors Ti(—imk n )~ l l 2 in front of each component {n^o-} provide conve- 
nient normalization. To connect bound Hamiltonian dynamics and scattering theory 
one should make the following very important step. Cut one part of CS off along CSOS 
(y = 0) and attach a semi-infinite waveguide (|5|) instead. Thus we introduce two scat- 
tering Hamiltonians 



-(h 2 /2m)d 2 + H'(y); ay>0, 
-(h 2 /2m)dl + H'(0)] ay<0. 



Let *& a (x,y, E) denote a scattering wavefunction which is a solution of the Schrodinger 
equation with the scattering Hamiltonian @ while it is also the solution of the ordinary 
bound Schrodinger equation (H) on the a— side C a (ay > 0). In the waveguide (ay < 0) 
any scattering wavefunction may be written in a form (0) where the two SOS states are 
no longer arbitrary but they are related through a generalized multichannel scattering 
operator T a (E) G C 

\-&- a } = f (E)\& a ) (9) 

since scattering wavefunction ty a (x,y, E) should satisfy proper boundary conditions on 
the er-side, ^ a (q G B a , E) = 0. 

Now one can easily give necessary condition for an energy E to be an eigenenergy of 
the original bound Hamiltonian H. Then an eigenfunction ty(x,y) should exist and the 
two scattering wavefunctions ty a (x,y, E) with the following SOS representation inside 
the waveguide 

VJx, y, E) = - 1 t={x\k- l/2 (E) \ e ia ^ E)y + e- laR(E)y f a (E)} |tf CT } (10) 
V —im 1 J 

and which match the eigenfunction outside the waveguide 

V(x,y) = ty sign{y )(x,y,E). 
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One should further require that eigenfunction ^(x, y) and its normal derivative d y ^(x, y) 
written in terms of ^ a ( x , y, E) should be continuous on CSOS (y = 0) which results in 
two simple equations (1 ± Tj-)|i?|} = — (1 ± Tj)|i?j_} which in turn are equivalent to the 
condition that is a fixed point of an operator T(E) = T^{E)T^{E) 

f(E)\^} = \^}. (11) 

Even stronger result can be proved [p3| : singularity of 1 — T(E) is also sufficient quanti- 
zation condition if an energy E is not equal to the threshold E' n for opening of some new 
mode n, and in this case the dimension of null space of E — H is the same as dimension 
of null space of 1 — T(E). 

The product of generalized scattering operators T(E) may be interpreted as CSOS- 
CSOS propagator or QPM although it is not a unitary but a compact operator and only 
the open-open parts of scattering operators T a (E) (consisting of matrix elements between 
open modes) are unitary matrices |23], |3(| . It is convenient to define wave operators which 
map any SOS state l^} to the corresponding scattering function outside the waveguide 
^ a {x, y, E) and vice versa 

(x,y\Q a {E)\^} = e(ay)* a (x,y,E), (12) 
mP a (E)\x,y) = 6(ay)*l(x,y,E*). (13) 

(0(y) is the Heaviside step function) where the dual SOS state {$*|, which is in general 
different from {^l, generates complex conjugated scattering wavefunction 

%(x,y,E*) = - 1 t={d* a \ \e iak{E)y + e- luk{E)y f a (E)\ K~^ 2 (E)\x}. 
y—ira L J 

The linear operators Q a (E) from C to TC and P a {E) from TitoC may be interpreted as 
the quantum CSOS-CS and CS-CSOS propagators, respectively. If E is an eigenenergy 
of the original bound Hamiltonian H and |t?|} is the associated fixed point of QPM 
T(E) then the corresponding wavefunction can be calculated by means of CSOS-CS 
propagator 

V(x,y) = (x,y\Q sign{y) \$ sign(y) } = (x,y\Q^} + faylQ^rf 

where := T^E)^^}. The position matrix elements of the energy-dependent quan- 
tum scattering propagators G a (E) = (E — H ff + iO)^ 1 will be written as G a (q, q', E) = 
(q\G a (E)\q'). Then we consider another, hybrid representation of scattering Green 
functions G a (y,y' , E), which are operator valued distributions acting over reduced SOS 
Hilbert space £ defined by 

{x\G a (y,y',E)\x'} = (x,y\G(E)\x',y'). 

The three newly defined quantum propagators can be expressed in terms of scattering 
resolvents |2]2], £3) 

f a {E) = —K 1 ' 2 {E)G a {^E)K 1 ' 2 {E)-l, 
m 

(x,y\Q*{E) = -^=0{ay){x\Cf a (y,O,^ 1/3 (^), 
P a {E)\x,y) = -±=0{ay)K 1 P(E!)£i l ,{0,y,Ef)\x}, 
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If one defines additional conditional CS-CS propagator (without crossing CSOS) Gq(E) 
as 

f (x,y\G^E)\x',y'); y>0,y'>0, 
(x,y\G (E)\x',y') = (x,y\G l (E)\x , ,y'); y<0,y'<0, (14) 

I 0; yy> < 0. 

then the usual {unconditional) energy-dependent quantum CS-CS propagator G(E) = 
(E — H)~ l can be decomposed in terms of four newly defined (CS/CSOS)-(CS/CSOS) 
propagators 

G(E) = G (E) + Y,Q°( E )( l - f -°( E )T°( E )Y l P-°( E ) 

+ Y^Q^^-T-^T^Y^-^P^E). (15) 

a 

All poles of the propagator G(E) — eigenenergies of H come from singularities of the 
form (1 — T(E))^ 1 so the quantization condition ([TT| ) is indeed justified. One is tempted 
to expand the factors like (1 — T(i?)) _1 into geometric series and give all the terms and 
factors a firm physical interpretation in terms of probability amplitudes for quantum SOS 
propagation [f[9| , p3[j . But there is serious difficulty since such sum would be manifestly 
divergent since the operators TL cr (£ , )T cr (i?) have in general also eigenvalues whose mag- 
nitude is (slightly) larger than one. Now we proceed to show that there exist another 
unique realization of quantum SOS propagators with unitary QPM and conditionally 
convergent expansion of decomposition formula. 



2.2 Reactance operator formulation and unitarization 

In this subsection we will unitarize the compact scattering operators T a (E) in three 
steps: (i) We will express the compact CSOS-CSOS propagators T a (E) as the Caley 
transformation of reactance operators [TBI], (ii) These reactance matrices will be made 
Hermitian by a simple transformation which preserves their open-open part (and thus 
preserving also the semiquantal and semiclassical limit) while it rotates the phases of 
components referring to close modes by 7r/4. (iii) Finally, we shall define unitary CSOS- 
CSOS propagators as the Caley transformation of Hermitian reactance operators. 

The solution of the Schrodinger equation (Q) is unique in any range of CS and at 
any energy when all the boundary conditions are known. Take any L 2 (5 ) function 
over CSOS (p(x) = {x\cp} and denote by ^ a (q,E) its extension as the solution of the 
Schrodinger equation (§) on the a— side C a of CS with the boundary conditions 

$ <T (x,0,E) = <p(x), $ a (qeB a ,E) = 0. (16) 

If on the other side of CS one assumes $ CT (q G C_ CT ) = this may be written as a linear 
relationship W a (E) from C to 7i 

\§)=W a {E)\ V }. 

For any value of energy E there exist nontrivial solutions of the Schrodinger equation 
on both sides of CS $> a (q,E) which have given values on the CSOS <f{x). The entire 
wavefunction Q(q,E) = $>^(q,E) + $±(q,E) is thus continuous on CSOS So by con- 
struction. But only for special values of E, eigenenergies, there exists such <p(x) that 
also the normal derivative of the wavefunction d y <&(x, y, E) is continuous at y = so 
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that &(q, E) becomes the eigenfunction — the solution of Schrodinger equation on entire 
CS C. The quantization condition d y ^(x,0,E) = d y <&i(x,Q, E) can be written as the 
singularity condition for the sum of operators over C 

(ME) + A^E))^} = (17) 

which are defined by 

{x\A a (E)\<p} = ad y $ a (x,0,E). (18) 

It is easy to see that the quantization conditions ( PD and are equivalent, that \<p} 
is related to |# CT }, A a (E) are related to T a (E), and W a (E) are related to Q a (E). From 
the scattering ansatz ( [T0| ) we see that one should put 

\^} = (i + f a (E)y 1 K l i 2 {E)\ v } (19) 

since then the scattering wavefunction ty a (q,E) becomes proportional to the wavefunc- 
tion $ CT (g, E) on C a . Differentiating fllCf ) with respect to y at y = one further obtains 

A a (E) =ik 1 l 2 (E)(l-f a (E))(l+f CT (E))- 1 k 1 l 2 (E) (20) 

according to definition (p!8]). Moreover, Q a (E)\ , d (T } = W a (E)\<p}, so 

W a (E) = Q„(E)(1 + Zm-'K^iE). (21) 

Relation ( p0|) calls for introduction of generalized (nonhermitian) multichannel reactance 
operators |18| 

R a {E) = k- 1/2 (E)A a (E)K- l/2 (E) (22) 

in terms of which one can write the generalized scattering operators as a Caley transfor- 
mation 

f a (E) = (1 + iR a (E))(l - iRv(E))- 1 (23) 

But for real energy E = E* the operators A a (E) are Hermitian. To see this one should 
transform surface integral over Sq to a volume integral over C a 

{<p\A a {E)\<p'} = o J dx$l(x,0,E*)d y K(xAE) = 
s 

= Jdq (d y K(q, E*)d y K(Q, E) + K(Q, E*)^(H'(y) - E)K(q, E)) . (24) 



Note that the second term of ( p4[) is also symmetric since H'(y) does not involve derivative 
d y . The operator K 1 / 2 ^) and therefore the reactance operators R a (E) are nonhermitian 
since the square roots of wavenumbers are not all real due to existence of closed channels. 
But there is an easy an unique way (up to trivial constant similarity transformation) of 
how to make operator K 1 ^ 2 ^) Hermitian without touching the components referring to 
open modes C a (E) which govern the physics in classically allowed regions of phase space. 
Let us define an operator 



2m 

*W = E \/ jr(E - E' n )sign(ReE - E' n )\n}{n\ (25) 



S 



which is a Hermitian for real energy E while it may be analytically continued for complex 
E for KeE ^ E' n . Its square root may be conveniently written in terms of simple piecewise 
constant unitary transformation K}I 2 (E) = u(E)K 1 ^ 2 (E), 

ReE>E' n ReE<E' n 

u(E) = \ n }i n \ + S v / isign(ImE - 0)\n}{n\. (26) 

n n 

We always consider the branch of square root with positive real part. We shall usually 
omit energy dependence of operator u since it is constant inside a given semiband E' n < 
KeE < E' n+l , ±lmE > 0. Now one can define Hermitian reactance operators in analogy 
with flU) 

K a (E) = t- 1/2 {E)A rT {E)t- l / 2 {E) = u-'R^u- 1 (27) 
and in analogy with (1231) one can define unitary operators 



%{E) = (1 + iK a {E)){l - i-R a {E))-\ (28) 

Unitary operators T^E) are no longer proper scattering operators of scattering systems 
(PD but they may be called unitarized scattering operators. The quantization condition 
for original Hamiltonian H (]17|) may be restated as the singularity condition for the sum 
of Hermitian reactance matrices 71(E) = 1Z^(E) + 7Z[(E) 

7Z(E)\p} = 0, (29) 

with ([I7|J27|) 

\p} = K}' 2 {E)\ v }. (30) 

Equivalently, this can be written as a fixed point (eigenvalue 1) condition for the product 
of unitarized scattering propagators T~(E) = T±(E)7^(E) 

f(E)\^} = \^} (31) 

since one can use the definition ( p8|) to derive a relation which connects ( p9|) and (|3l|) 

1 " f{E) = 1(1 + 71(^)^(^(1 + 7}(E)) 

and the two types of stationary SOS states are related by 

\p} = \*h} + \1>l}, m = ^ + %{E)Y 1 \p} (32) 
where we have defined \ipi) := T^(E)\tp^} so that general relation analogous to @ holds 

\^ a } = %(E)\^ a }. 
The fixed points of compact and unitary QPM are related by the formula 

K} =[iL- \{U - 0(1 + Zr(E))] |^}. 

which follows from eqs. ( |19| , ^0| ^). 
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2.3 Unitarized decomposition formula and interpretation 

Now we show that we can interpret the operator T(E) as CSOS-CSOS propagator or 
QPM which has all desirable properties. One can define also the unitarized versions of 
other three conditional CS/CSOS-CS/CSOS propagators 

Q a (E) = Q a (E) [u-Ku-u-^il + f^E)) 

VAE) = [u - ±(1 + -iT 1 )]^) (33) 

Go(E) = G (E)+ l -^2Q a (E) [l-u*-l(u-u-i)(l + f a (E))(u-u~ l )]P a (E) 

a 

and show by simple algebraic manipulation that the unitarized decomposition formula 
for energy-dependent quantum CS-CS propagator follows from (p~5|) 

G(E) = g (E) + Y.Q<r{E){l-f_ a {E)%{E))- 1 V^{E) 

(7 

+ YsQ^^-T-^ME^f^E^E). (34) 



In definitions ([33]) we have multiplied by simple linear combinations of unitary operators 
so we have not introduced any singularities for real E. All poles of the resolvent G(E) 
again come from singularities of (1 — T(.E)) _1 in accordance with (|31|). If j^t} is a fixed 
point of T~(E) for such an eigenenergy E then the corresponding wavefunction is given 
by 

V(x,y) = (x,y\Q S i g n(y)\^si g n(y)} = (x, y\ Q T (£) |^ T } + {x, y\ 

where '■— 7f (E)\i/)^}. One can expand the factors (1 — T(E))^ 1 in unitary decom- 
position formula in a geometric series 

oo 

<?(£) = (35) 

n=0 

with 

$2i+i{E) = ^&(B)(f-„W? a (E)) l f-„(E!)P (T (E), (36) 
G 2l (E) = Y.^{E)(f- (T {E%{E)) l V^{E), 1 = 0,1,2... (37) 

a 

This decomposition formula may be given the following firm physical interpretation. 
The quantum probability amplitude (q\G(E)\q') to propagate from point q' to point 
q in CS at energy E may be written as a sum of conditional probability amplitudes 
(q\Qi(E)\q r ) to propagate from q' to q at energy E and cross CSOS exactly I times (f55|). 
If points q and q' lie on the same/opposite side of CS with respect to CSOS then each 
continuous orbit must cross CSOS even/odd number of times and thus all probability 
amplitudes (q\Qi(E)\q') for odd/even I are zero. Each probability amplitude (q\Qi(E)\q') 
may be further decomposed as a sum of products of probability amplitudes by inserting 
an identity J c?a;|a;}{a;| between each pair of factors in ( |36| , |37D . The elementary prop- 
agators should be interpreted as follows: (q\Q a (E)\x'}, {x\V a {E)\q'), {x\f a (E)\x'} are 
conditional quantum probability amplitudes for a system to propagate at fixed energy 
E through a— side of CS C a from point q' in CS or point x' on CSOS to point q in CS 
or point x on CSOS and without crossing CSOS in between. 
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Each orbit of a bound system which crosses CSOS must cross it again with probability 
one, consistently with unitarity of CSOS-CSOS propagators 

J dx\{x\%(E)\x'}\ 2 = l. 
2.4 Quantum Poincare evolution 

Let us choose normalized initial SOS state |0, f}. {x\0, f} should be interpreted as 
quantum probability amplitude that system's orbit initially crosses CSOS at x from 
below. Like for classical Poincare evolution one should also specify the value of energy E 
besides CSOS coordinates x to completely determine system's dynamics. Then quantum 
Poincare evolution is a simple iteration of unitary QPM T(E) 

|n + l,T} = t(£)KT}. (38) 

When the system crosses CSOS n-th time from below, {x\n, f} = {x\T n (E)\0, f } is a 
probability amplitude that it crosses CSOS at x if it was initially in a state |0, |}. Ex- 
istence of a stationary state implies existence of a fixed point of QPM T(E) which does 
not evolve, not even its phase. 

Resume: I have shown that the /-dim autonomous bound Hamiltonian quantum dy- 
namics governed by a continuous group of unitary propagators {exp(—itH/h); t G 71} is 
equivalent to a continuous family (labeled by energy) of (/ — l)-dim discrete quantum sys- 
tems governed by discrete groups of unitary propagators {T n (E); n = 0, ±1, ±2 . . .} just 
like the classical autonomous bound Hamiltonian dynamics with / freedoms is equivalent 
to a continuous family (labeled by energy) of discrete area preserving classical Poincare 
mappings with (/ — 1) freedoms. 



3 Phase space SOS representation of quantum SOS 
states 

For some purposes one would like to have as close correspondence with classical Poincare 
dynamics as possible. The classical system can be efficiently described in terms of prob- 
ability distribution function over the 2(/ — l)-dim SOS with coordinates z = (x,p x ). 
Classical Poincare mapping te on points z or distribution functions f(z) can be written 

as 

z n +i = r E (z n ), (39) 

fn+l(T E (z)) = fn(z). (40) 

Classical invariant distributions, fixed points of te, are characteristic functions over 
disjoint classical invariant components of SOS which can be either periodic orbits, or 
(/ — l)-dim invariant tori, or 2(/ — l)-dim chaotic regions. From the full space quantum 
mechanics we adopt coherent state representation for the quantum SOS distribution of 
an SOS state \ip} and call it SOS Husimi distribution 

^ a (z) = \{z,am\ {x'\x, Px ,a} = (2nha)^e W ^- { -^^ + ^^y (41) 
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which is not unique but depends on the deformation a of the minimal coherent states 
| z, a}. In completely similar way as for ordinary continuous time ( / — 1) — dim quantum 
system one could also study the familiar Wigner ]33| distribution but for the purposes of 



this paper I prefer Husimi distribution because it is positive unlike Wigner. According 
to most general and strongest correspondence principle, the so called principle of uni- 
form semiclassical condensation || |29], in the semiclassical limit h — > phase space 
distributions of quantum eigenstates as well as quantum states evolving for sufficiently 
long time should condense on classically invariant components. So SOS Husimi distribu- 
tion ^f^(z) of a typical eigenstate, where \ip} is the corresponding fixed point of unitary 
QPM, should appear very much like classical plot of chaotic region or some regular in- 
variant torus. In the next section we present numerous numerical results concerning the 
demonstration of this principle in a generic nonlinear dynamical system obtained by the 
quantum SOS method. 

For quantitative analysis we define localization area (volume) of a given normalized 
SOS state \ip} through its information entropy 

Ap = c/_iexp (- j dz^ h a (z)M h a (z)) (42) 



The localization area A$ measures the effective area of SOS occupied by a state \ip} i.e. 
the area where the SOS Husimi distribution is significant. c/_i is some dimensionless 
normalization constant which may be determined by the requirement that Husimi func- 
tion generated by truly Gaussian random wavefunction, the so called Gaussian Random 
Husimi distribution (GRHD), should have localization volume equal to the volume of 
entire classically allowed region of phase space, yielding numerically in 2-dim c\ ~ 1.538. 
In order to get more information about the structure of SOS states one could use the 
concept of generalized entropies to define corresponding generalized localization areas 
A^is) for s > 



V) = °/-i(«) /<**(*£(*)) 



(43) 



with the usual localization the limit = A^(s — > +0). Again we may use 

GRHD to determine normalization constants c^_i(s). Here are some numerical values 
ci(l) « 2.02, ci(2) « 2.49, ci(3) « 2.94. 



4 Application of quantum SOS method to quantum 
chaos 

In this section I will present the results of extensive numerical application of quantum 
SOS method to a generic 2-dim bound autonomous Hamiltonian system, namely the 
semi-separable oscillator. Due to efficiency of quantum SOS method which effectively 
reduces the labor by one degree of freedom and special geometric structure of our sys- 
tem we were able to reach extremely deep semiclassical regime, namely we were able 
to calculate thousands of consecutive levels with sequential quantum number around 
twenty million. Thus we are able to test various conjectures of quantum chaos about the 
structure of eigenstates and statistical properties of spectra. 
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4.1 The system and practical quantization technique 

The most practical quantization condition of quantum SOS method is the singularity 
condition for the sum of Hermitian reactance operators (p9|) 

det R(E) = 0, (44) 

which are for numerical calculations represented by finitely dimensional matrices R(-E) = 
R-f(-E') +Rj(-E) in a truncated (N = N Q + N c )-dim basis of N Q = dim C a (E) open modes 
and sufficiently large number N c of first closed modes such that the zeros of flH]) converge. 
Let ty an (x, y, E) denote the two unique sets of solution of Schrodinger equation (|2|) on 
either a— side C a and with boundary conditions ^ an (x, 0, E) = {x\n}, ^ an (q G B a , E) = 
0. In other words, \^/ an (E)) = W a (E)\n}. For real energy E the matrix elements of 
Hermitian reactance matrices can be written as 

R anl (E) = {n\K(E)\l} = [ dxV* an {x,y,E)d y y al {x,y,E)\ ay=+0 . (45) 

■ /\k n (E)ki(E)\ J 



If the system possesses a time-reversal symmetry then the Hermitian reactance matrices 
R CT are real (and symmetric) due to reality of wavefunctions ^ an (q, E) and therefore 
their Caley transforms, unitary conditional CSOS-CSOS propagators have symmetric 
matrices {n\%(E)\l} = {l\%{E)\n} . 

Reactance matrices can be most easily calculated for the so-called semiseparable 
systems, that is, for systems which are separable (in (x, y) coordinates) on both sides of 
CS C a but they have possible discontinuity on CSOS so that they are not separable on 
the whole CS C 



One such semiseparable system which turned out to be very convenient for numerical 
work is the so-called (/ = 2) -dim semiseparable oscillator (SSO) with the Hamiltonian 

H(x,y,p XJ p y ) = \p 2 x + \p 2 y + l(x+ |sign(y)a) 2 , -b l < y < 6 T (46) 

The potential is harmonic in x— direction while it is flat with perfect hard walls at y = 
— in y— direction. The classical dynamics of SSO was also extensively studied and 
it was found that the system exhibits all the features of generic nonlinear softly chaotic 
2-dim autonomous Hamiltonian systems. For limiting cases a = (single box limit) and 
a = oo (two box limit) the system is integrable, while for most other values of parameters 
the system has mixed classical dynamics with regular and chaotic regions coexisting in 
phase space and on SOS (x,p x ;y = 0), for some values of parameters the system is 
even fully chaotic - ergodic (see e.g. figure 5c). The quantized SSO described by 
Hamilton operator H = H(x, y, —ihd x , —ihd y ) has a 1-dim scaling symmetry (a, b a , E) — > 
(Aa, Xb a , X 2 H, X 2 E). The reduced Hamiltonian is just a simple 1-dim harmonic oscillator 
H' = —\% 2 d 2 + \x 2 with real SOS eigenmodes 

{x\n} = {n\x} = (V^h2 n n\)- 1/2 exp(-x 2 /2fi)H n (x/Vh) 
and threshold energies E' n = (n + |)fr determining the wavenumbers 



k n (E) = hr l J2E - {2n + l)h, n = 0,1,2 



3 The ergodicity in such special cases has not been rigorously proved but it can be shown numerically 
that the total volume of regular components can be made unmeasurably small. For example, for any 
< a < \/8i? the system becomes ergodic in the limit + — > oo. 
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Due to separability the solutions of Schrodinger equation on each side C a are composed of 
products {x + \aa\n} sm(k n (E)(b a — cry))/ sm.{k n (E)b a ) but because of defect on CSOS 
at y = one should use unitary shift operator O = exp(iap x /2h) with matrix elements 

On/ = J dx{n + ^a\x}{x\l} 

to generate SOS induced solutions 

where we write = O n i, 0[ n i = 0/ n . So, using (|45"1), real symmetric reactance 

matrices for SSO read 

R anl (E) = -\k n {E)ki{E)\- 1 / 2 Y, -nM E )^t{k j {E)b a )O a ^ (48) 

3 

The matrix elements of shift operator can be numerically calculated via stable symmetric 
recursion 

I / .2 
O n , = (-l) n O ,„ = ^=exp 



O r 




n,l-l- 



The number of open modes of SSO is N a = round(E/h). In order to determine minimal 
number of closed modes N c , such that the results are expected to converge one can use 
semiclassical arguments, namely, SOS phase space supports of coherent state represen- 
tation of basic SOS states \n}, < n < N a + N c — 1 should cover the supports of first 
iV shifted states 0\n}, < n < N Q - 1 

N c = —^= H \N . 

\V2E 2E 



It is very important to stress that the shift matrix O n i and therefore also the reactance 
matrices tl an i(E) are effectively banded, that is, their elements are decreasing expo- 
nentially fast when the distance from diagonal \n — l\ becomes larger than the effective 
bandwidth. One can again derive semiclassical formula for the effective bandwidths us- 
ing overlap condition for the coherent state representation of the SOS states \n} and 
0\l} 

2a 

bandwidth (R CT (E)) = 2bandwidth(0) w ^=N Q (49) 



'2E 

Note that the function f[E) = detR(i?) has singularities (poles) at the points E where 
for some n, k n (E)b a is a multiple of ir. But between the two successive poles f(E) is 
smooth (even analytic) real function of real energy E. I have devised an algorithm for 
calculation of almost all levels — zeros of f{E) within a given interval [Ei, Ef] which 
needs to evaluate f(E) which takes O (bandwidth 2 N) FPO only about 25 times per 
mean level spacing while it typically misses less than 0.5% of all levels. The control over 
missed levels is in general a very difficult problem. The number of all energy levels below 
a given energy E, M{E) can be estimated by means of the Thomas-Fermi rule 

Af(E) w Af TF (E) = b -l±^k(2E)V 2 . (50) 

37m 
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But this formula is generally not very helpful even if next semiclassical corrections are 
negligible since the fluctuation of the number of levels in an interval [Ei,Ef] is propor- 



tional to \jN(Ef) — N(Ei) except in the extreme case of fully chaotic systems where 
the spectra are much stiffer and the fluctuation is proportional to \og[J\f(Ef) — Af(Ei)] 
so that Thomas- Fermi rule can be used to detect even single missing level ]7|. |l[ . 

We have chosen the following values of parameters a = 0.03, bi = 5.0, = 10, E — 
0.5, while for quantal calculations we take the energy to be in a narrow interval around 
E = 0.5. Careful examination of classical dynamics (for Poincare SOS plot see figure 4w) 
showed that there is only one dominating chaotic component of phase space of relative 
volume p2 = 0.709 ± 0.001 (which is not equal to its relative SOS area) while regular 
component together with other very small chaotic components have total relative volume 
pi = 1 — p2 = 0.291 ± 0.001. We have calculated two stretches of consecutive energy 
levels and corresponding eigenstates: (case I) 14 231 levels in the interval 0.35 < E < 0.65 
for h = 0.01 with sequential quantum number according to ( [50"D equal to M ~ 16 000, 
(case II) and 13 445 levels in the interval 0.49985 < E < 0.500105 for % = 0.0003 with 
sequential quantum number Af ~ 17684 000. 

Large square root number fluctuations prevent to determine the number of missed 
levels by using Thomas-Fermi rule (although higher order semiclassical corrections are 
negligible in this regime). One can compare the number of levels M{E) with the number 
of levels Afo(E) or A r 00 (-E') for the two nearby integrable - separable cases (with the 
same b a but with a = (single box limit) or a —>■ oo (two box limit), respectively) since 
the leading order semiclassics (Thomas-Fermi rule) does not depend upon the defect a. 
Afo(E) and Noo(E) can be easily calculated numerically and large scale fluctuations of 
M{E) — J\fo,oo{E) turn out to be much smaller than the fluctuations of J\f(E) —Af TF (E) 
suggesting that we have missed less than 20 levels out of 14 231 at % = 0.01 and 40 - 80 
levels out of 13 445 at h = 0.0003. Note that in the first case (h = 0.01) there was much 
less almost degenerate pairs of levels (and therefore less missed levels) due to the level 
repulsion. 

For each zero of equation (fBJ), eigenenergy E, one can determine the components 
pn = { n \p} of SOS representation \p) of the corresponding eigenstate by solving the 
homogeneous equation 

£R n/ (£)# = 0. 
i 

The corresponding eigenfunction ty(x,y) can be written as 

(a;, y) = £ - f £!==^ gaiy)l (x, y, E) (51) 

\H E )\ 



where the wavefunctions y, E) are given by (|4"7"D. For the SOS Husimi distribution 

of (eigen)states of SSO it seems natural choice to take coherent states with a = 1 since 
basic SOS state |0} is then just a coherent state located at the origin z = (0, 0). In polar 



coordinates, r = Jx 2 + p 2 ,, (p = arctanQo^/a;) the corresponding SOS Husimi function 



2nh 



u ■ -in<j> 



2 



e , (52) 
V2h 



n Vn! 

can be efficiently calculated by means of Fast Fourier Transformation where ip n = {^I^t) 
are the components of the corresponding fixed point of the unitary QPM = |(1 — 
in^)\p} (see eqs. (|3|j2|), _^ 

Ipn = Pn - i E ^nl{E)pi- 

l 
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4.2 Classification of eigenstates and Berry-Robnik level spacing 
distribution 



According to principle of uniform semiclassical condensation the phase space distribu- 
tions of eigenstates (such as Husimi) should uniformly condense on the classical invariant 
components of phase space, which may be either regular — tori, or chaotic, when h — > 0. 
This condensation should be understood in a weak sense, i.e. quantum phase space 
distribution smoothed over many Planck's cells should approach characteristic function 
of a classical invariant component. Thus any quantum state, if % is sufficiently small, 
can be classified either as regular or chaotic, if it is associated with classically regular or 
chaotic component, respectively. 

Berry and Robnik Q also assumed that the energy levels of states associated with 
different disjoint classical components cannot be statistically correlated, so that the en- 
tire energy spectrum is a superposition of statistically uncorrelated level subsequences 
associated with different classical invariant components. The most characteristic is the 
so-called level spacing distribution P(S), where P(S)dS is a probability that a randomly 
chosen spacing between two adjacent energy levels lies between S — dS/2 and S + dS/2. 
All regular levels may be merged together giving a totally uncorrelated sequence with 
the so-called Poissonian statistics with -Ppoisson(S') = e~ s while each chaotic subsequence 
is statistically equivalent to the spectrum of a fully chaotic system and therefore also 
II to the spectra of infinitely dimensional Gaussian orthogonal /unitary random matri- 
ces (GOE/GUE) provided that chaotic states are mainly delocalized - extended over 
the whole chaotic component. The gap distribution E(S) = fg 2 da(o~ — S)P(a) factor- 
izes upon statistically independent superposition of spectra, so if one assumes only one 
practically dominating chaotic component with relative volume p2 and regular and tiny 
chaotic components with total relative volume p\ = 1 — p 2 the ultimate semiclassical 
Berry-Robnik formula reads 

E™(S) = E p —( Pl S)E GO *(p 2 S), P™(E) = -^E^(S). (53) 

This two component Berry-Robnik formula will apply also to SSO with a = 0.03, 6-f = 
5, 6j = 10, E — 0.5, where we have indeed only one large dominating chaotic region 
(see figure 4w). The Berry-Robnik distribution does not exhibit level repulsion, since 
P^ R (0) = 1 — p\ 7^ 0. On the other hand there has been a vast amount of phenomeno- 
logical evidence [^8j in favour of the so called fractional power law level repulsion which 
is globally very well described by the Brody || distribution 

Pf{S) = a^exp(-65 /3+1 ), a = {(3 + 1)6, b = [T(l + ((3 + l)" 1 )]^ 1 (54) 

which is characterized by the noninteger exponent (3, P(S — > 0) oc S 13 . Numerical spectra 
which contain even up to several ten thousands energy levels of quantum Hamiltonian 
systems with mixed classical dynamics typically still exhibit the phenomenon of fractional 
level repulsion, with statistically significant global fit by the Brody distribution. In such 
cases there was a persisting puzzle as for how the level spacing distribution converges 
to the semiclassical Berry-Robnik distribution as one increases the sequential quantum 
number or decreases the value of effective h. However, recently we have succeeded 
to demonstrate the ultimate semiclassical Berry-Robnik level spacing distribution in a 
rather abstract 1-dim time-dependent dynamical system, namely the standard map on 
a torus, and showed (smooth) transition from Brody-like to Berry-Robnik distribution 
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as h decreases |27|, |28| (see also |2T|), and more recently the first such demonstration in 



a generic 2-dim autonomous conservative system is provided by the SSO |24 . 

Let us estimate the maximal (critical) value of the effective Planck's constant h max of 
the far semiclassical regime where Berry-Robnik approach is expected to be valid. There 
are two conditions to be satisfied: 

• States which live on classically disjoint invariant components should have small 
overlap to provide statistical independence of partial subspectra. Husimi phase 
space distributions of quantum (chaotic or regular) eigenstates typically decay 
like Gaussian into classically forbidden neighbouring (regular or chaotic) invariant 
region with an effective penetration depth equal to \/h (see e.g. [0). It seems 
reasonable to require that this quantum resolution scale should be much smaller 
than dimensions of classical chaotic region, say at least 10 times smaller than the 
radius of the largest ball which lies entirely in the chaotic component of SOS, 
\/h ~ 0.1r c . For SSO with a = 0.03, — 5, — 10 we have r c ~ 0.16 (see figure 
4w), so h ~ 0.0003. 

• Chaotic states should be delocalized - extended over the whole classical chaotic 
region of phase space in order to justify usage of maximal entropy ensembles of ran- 
dom matrices (GOE/GUE) to model chaotic subspectra. If ^ is a relative area of 
chaotic component of SOS then SOS Husimi distributions of last /^^o open SOS 
eigenmodes \n} (thin circular rings in case of SSO) approximately support the 
chaotic region (which has the shape of a ring in case of SSO (see figure 4w)). Quite 
generally J25[] QPM is represented by a banded matrix (having approximately inde- 
pendent /i2-?Vo-dim chaotic block) with a minimal bandwidth b « Aa^No/ y/2jI^E 
in case of SSO. Using the theory of localization of eigenvectors in finite banded 
random matrices [12| one can write the condition for, say 90% derealization of 



chaotic states of SSO 1.4-^ ~ i-'o 9 ' giving 

< „2 



h ~ a 2 . (55) 



For SSO with a = 0.03 we have the condition h ~ 0.0009. 
The far semiclassical regime sets in if 

h ~ ftmax = min(0.01rc, a 2 ) 

and the case II has been chosen to optimally meet this condition. 

Unique classification of states into regular and chaotic class is a necessary condition 
for the validity of Berry-Robnik formula and it can be performed quantitatively as fol- 
lows. Let us define a probability distribution of localization areas V(A) where V(A)dA is 
a probability that a randomly chosen eigenstate have SOS localization area (f42|) between 
A — dA/2 and A + dA/2. For two component case with a single dominating chaotic region 

V(A) is expected to be bimodal, one sharp peak of width proportional to O ^h^~^j close 

to A = correspond to regular states while the other, wider peak located at SOS area of 
chaotic component correspond to chaotic states. The width of the second peak is a more 
complicated and yet unsolved function of h and geometry. Find the area A m[n between 
the two peaks where V(A) takes its minimum. A quantum state with SOS localization 
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area A is said to be regular if A < A min and chaotic if A > A min . The integrated proba- 
bility for a randomly chosen state to be regular should match with the relative volume 
of classical regular component 



pf= J dAV(A) (56) 
o 

This was nicely confirmed for the stretch of 14 thousand eigenstates of far semiclassical 
case II (h = 0.0003), giving p^ = 0.293 ± 0.004 in excellent agreement with the classical 
value pi = 0.291 (see figure la), whereas for the near semiclassical case I (h — 0.01) 
the localization area distribution V(A) was still unimodal, indicating that eigenstates 
cannot be clearly and uniquely classified as regular or chaotic (figure 2a). 

The level spacing statistics behaves in complete agreement with the results of this 
classification. For case II one obtains significant fit (x 2 = 12150) with Berry-Robnik 
distribution (|53"D (figure 3b). I have also separated the entire spectrum into regular 
and irregular part according to classification in terms of localization area statistics and 
studied level spacing distribution for each part separately. The statistics of the regular 
part which contains 3791 levels, is indeed very close to Poissonian (figure 3c). The 
statistically significant fit with Berry-Robnik distribution gives p\ = 0.86, x 2 — 1500 
while fit with Brody distribution gives f3 = 0.006, x 2 — 2300. The statistics of the 
irregular part, which contains 9652 levels, is slightly further from but still close to GOE 
(figure 3d). The Berry-Robnik fit is only slightly nonsignificant, pi = 0.04, x 2 — 15000, 
while Brody fit is worse, giving P = 0.83, x 2 = 23000. These results (figure 3) clearly 
confirm Berry-Robnik picture (which is claimed to be an asymptotically — as h — > - 
exact theory [f26j [27], |28| ) although we still see small but still significant deviations from 



the expected statistics of partial spectra, because of small but still existing correlations 
between regular and irregular levels due to (small) overlap of corresponding eigenstates, 
and because of localization of some irregular states on small subregions of chaotic orbit 
where classical dynamics is almost trapped, such as near chaos border (see also figures 
4i,4t)). On the other hand, energy level spacing distribution of the near semiclassical 
case I still exhibits power law level repulsion with statistically significant fit by Brody 
distribution (|5J) with (3 = 0.142, X 2 = 5320 (figure 3a). 

There is another useful quantitative measure of a given SOS state namely, 
quantum-classical overlap: the overlap between SOS Husimi distribution and the classical 
chaotic component Sc of SOS 



Bjf, = / dz^ a (z). (57) 
Js c 

The state is irregular if is close to 1 and regular if is significantly less than 
1. Similar classification has been recently performed in a Robnik billiard Again 
we define the probability distribution of quantum-classical overlaps V(B): V(B)dB is a 
probability that the quantum-classical overlap of a randomly chosen state lies between 
B — dB /2 and B + dB/2. We expect and confirm (figures 1 and 2) that this distribution 
has qualitatively the same properties (bimodality) with even sharper peaks as localization 
area distribution V(A). V(B) is bimodal even for the case I and thus provides a measure 
for classification of states for near semiclassical regime where localization area statistics 
fails. But for far semiclassical regime (case II) the V(A) has lower minimum than V{B) 
and thus provide clearer classification of eigenstates. The scatter diagram (figures lf,2f) 
A^ versus B^ is very interesting and helps us to discover states with some special or 
exotic geometry, e.g. the regular states which live on small regular islands inside chaotic 
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region (see figure 4(g,h,r,s,w)). These states have small localization area A whereas their 
quantum-classical overlap B is significantly different from zero while for all other regular 
states living in the large central island B is extremely close to zero. 

4.3 Gallery of eigenstates 

For the systems of two freedoms (/ = 2) SOS Husimi distributions provide certainly the 
most elegant and efficient way of graphical presentation of quantum eigenstates. In figure 
4 I present the SOS Husimi distributions of eleven eigenstates of far semiclassical case II. 
The first six of them are typical consecutive states within the qualitative conclusions of 
previous subsection: the two regular states are strongly localized over corresponding tori, 
while three out of four chaotic states are very much uniform over the chaotic component 
with lots of (GRHD-like) microscopic structure which is responsible for the normalization 
factor Ci = 1.54 in the definition of localization area (]42"1). The fourth chaotic state 
(figures 4(e,p)), which is rather exceptional in this regime, is mainly localized only in 
the outer part of chaotic component due to circular partial classical barrier at r « 0.81 
0. The remaining five states are the specially picked states which have geometrically 
extremal properties: (i) the regular states with smallest area and smallest quantum- 
classical overlap on the small regular islands inside chaotic component (figures 4(g,h,r,s)), 
(ii) the chaotic state which is strongly localized around the chaos border where the 
classical dynamics is almost trapped (figures 4(i,t)), (iii) the mixed state which is almost 
uniformly localized over chaotic component and a torus which is geometrically close 
to chaotic component (figures 4(j,u)), and (iv) a chaotic state which lives on a small 
broken separatrix (figures 4(k,v)). We always give two graphical representations of SOS 
Husimi distributions: equidistant contour plot with ten contours separated by 1/10 of 
the maximal value of SOS Husimi distribution starting with 1/20 of the maximum, and 
logarithmic five contour plot where neighbouring contours are by a constant factor of 
h~ ' 2 apart (the values of SOS Husimi distribution at the first and last contour differ by a 
factor of 1 /h). The two presentations are complementary: equidistant plot shows only the 
most important features of an eigenstate whereas logarithmic plot slightly obscures the 
most important features and shows also the important details of SOS Husimi distribution 
such as extension over classical invariant components or distribution of its zeros. 

I have also calculated a sequence of 16 consecutive eigenstates in a classically fully 
chaotic — practically ergodic regime a = 0.25, b-\ = 4, &^ = 11, % = 0.0003,0.49999 < 
E < 0.4999903 (case III). The sequential quantum number is here the same as for case II, 
M ph 17 684 000. SOS Husimi functions of all 16 states were very uniformly extended over 
the whole classically allowed region of SOS which confirm one of the basic conjectures 
of quantum chaos (see figure 5) ||. 

In figure 6 I show the CS wavefunctions of two typical eigenstates of case II, a regular 
and a chaotic, which were calculated according to formula (|51|). Despite extremely high 
sequential quantum numbers the chaotic wavefunctions are globally scarred with many 
(stretches of) classical orbits, whereas on smaller scales they appear much more Gaussian 
random. On the other hand, wavefunctions of fully chaotic case III appear uniformly 
random (without any structure) even on the largest scale (an example is given in figure 
7). 

In order to test numerical method and to search for possible large scale nonuniformi- 
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ties in wavef unctions I have denned also a contrast of a wavefunction ty(x,y) by 



* fd x fXdy(^ + e -^)Mx,yW [ } 

Again, 9{y) is the Heaviside step function. The contrast Cy lies between —1 and +1 
and measures the difference between quantum and classical probability that the system 
lies above or below CSOS. I have also studied the probability distribution of contrasts 
V(C) for cases I and II (figures lc,2c): V(C)dC is a probability that randomly chosen 
eigenstate has a contrast between C — dC/2 and C + dC/2. Contrast may be used to 
detect scars of periodic orbits. I found that large majority of eigenstates has a contrast 
very close to zero. But contrast may be significantly different from zero for some chaotic 
states, since contrast Cy is strongly correlated with localization area Ay and especially 
with quantum-classical overlap By (see scatter diagrams l(d,e),2(d,e)). These states 
must be scarred by classical periodic orbits since these are typically not balanced with 
respect to CSOS. It has been numerically checked and found that by far most frequent 
are the so called bouncing-ball-like scars, which can have the largest contrast (up to 
1), and which are associated with a continuous family of neutrally stable orbits with 
y = const. I show one example for case II in figure 8. 



4.4 Distribution of zeros of SOS Husimi distribution 

According to Bargmann the 2-dim Husimi distributions of quantum states may gener- 
ally be written as \ f(x+ip x )\ 2 exp(— (x 2 +p 2 )/2h) where / is a complex analytic function. 
This is indeed the case for SSO fl52"|). The zeros of SOS Husimi distributions are thus sim- 
ple points in SOS and their distribution can be associated with the dynamical properties 
of an underlying system |T7|. The zeros of SOS Husimi distributions of chaotic states 



are expected and confirmed to be uniformly distributed over the chaotic component of 
SOS (figure 9, especially if the system is fully chaotic (figure 10) P7|). SOS Husimi 
distribution of a chaotic eigenstate in a generic mixed system have also zeros on the 
regular component where the zeros condense on the regular invariant curves — tori (see 
figure 9). The large number of the zeros of regular states in a mixed system lies on some 
curves which are not classically invariant and which can even extend to chaotic region 
(generalization of anti-Stokes lines from ||17|| ) while substantial number of zeros are also 
uniformly distributed over chaotic region and along regular invariant curves. But zeros 
for a regular state typically strongly avoid the classically invariant region of strongest 
localization of Husimi distribution (figure 9(a,d,g,h,j,k)). Even chaotic states of a fully 
chaotic system can have (1-dim-like) clusters of zeros lying outside classically allowed 
region of phase space (see figure 10). 

So far there are no analytical or numerical results about other statistical measures 
of zeros of Husimi distribution of a chaotic state inside chaotic region, such as e.g. 
nearest zero spacing distribution Z(S): Z(S)dS is a probability that a distance between 
a randomly chosen zero and its nearest neighbour lies between S — dS/2 and S + dS/2. 
I have found numerically: (i) that zeros of SOS Husimi distribution which live in the 
chaotic region for / = 2 feel a cubic repulsion Z(S — > 0) oc Z 3 (figures 11,12), and (ii) 
the nearest zero spacing distribution for chaotic regions of chaotic states can be very 
well modeled by GRHD (the same as used to define the normalization Cf of localization 
areas) (figures 11,12), except for large spacings S it seems (see figures lib, 12b) that 
Z(S) behaves like hiZ(S) oc — Z 4 while GRHD model suggests a Gaussian behaviour 
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Z(S) oc — Z 2 . Further numerical and if possible analytical work is required to clarify 
and explain these observations. 



4.5 Demonstration of Quantum Poincare evolution 

The physical motivation for an explicit study of quantum Poincare time evolution in 
SSO is a quantum-classical correspondence: How long can a quantum evolution follow 
classical evolution? In order to explore this question is some detail we study a quantum 
Poincare evolution of initial wave packet — coherent state \z} = \z,a = 1}, z = (x,p x ) 
at fixed energy E in SOS Husimi representation 

h n (z f ,z t ,E) = \{z f \f n (E)\ Zl }\ 2 . (59) 

Classically, we should take the same initial wave packet SOS phase space distribution 
|{z|zj}| 2 , then classically evolve it fl4"0|), and finally take the classical probability that 
the system finds itself described by the final wave packet distribution |{z|z/}| 2 

f n (z f , Zi ,E) = J dz\{z f \z}\ 2 \{r^z)\ Zi }\ 2 . (60) 

Such smoothed classical evolution of a wave packet f n ( Z f, Zi ,E), which is sometimes 



called coarse-grained classical dynamics [33], is a purely classical (not semiclassical) ob- 



ject which is expected to be most faithfully followed by the quantum SOS phase-space 
propagator h n (zf, Zi ,E). I plot these quantum and classical SOS phase space distribu- 
tions for initial wave packet located somewhere in chaotic region for far semiclassical 
case II (figure 13) and near semiclassical case I (figure 14). In agreement with the results 
of stationary quantum mechanics we find, that in case II quantum dynamics explores 
the whole classically accessible part of SOS (the chaotic component) while in case I 
it remains localized on much smaller subregions of chaotic component. There exists a 
kind of break iteration ribreak up to which quantum dynamics faithfully follows classical 
dynamics |L0|. We can characterize this quantitatively by means of localization areas 



A^,A^ of quantum and classical SOS distributions, h n , f n , respectively. (Note that for 
calculation of localization area of classical SOS distribution one should take Cf — 1 in 
formula P2|).) So, for n > n break , A^ become significantly smaller than A^. But A^ may 
still be increasing for n > n^ rea ^ up to some n satU r where it saturates and then fluctu- 
ates around some average value A^. The far semiclassical regime, where the ultimate 
semiclassical formulas of quantum chaos (e.g. for level spacing distribution [ffl, for delta 
statistics |32| , or for the statistics of matrix elements |20|, ^TJ) are expected to hold, can 



be defined by the condition that arbitrary initial wave packet should explore the whole 
classically accessible region of phase space, i.e. 



AL = A 



This is equivalent to the condition that eigenvectors of QPM T(E) should not be localized 
inside classically invariant chaotic components of SOS (|55|) . In the opposite case QPM 
evolution of initial wave packet takes place only in the more or less small subspace 
spanned by these eigenvectors of QPM which have significant overlap with initial state. 
This is due to quantum localization |9], [L2| which is a consequence of handedness of QPM. 
Sufficient condition for the far semiclassical regime is that classical SOS distribution 
should reach an equilibrium within nb re ak iterations, A C J: rs A^L. 

j- ' "-break l -* J 
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To demonstrate these phenomena most clearly I have studied (in figure 15) the SOS 
dynamics of SSO in the so called diffusively ergodic regime of small a <C 1 and very 
large 6j + 6j_ 3> 1 where classical dynamics is ergodic although it explores the accessible 
phase space rather slowly, with r 2 = x 2 + p 2 being an approximate second integral of 
motion. For a = 0.03, 6-f = 500, foj = 1000 the classical dynamics equilibrate after several 
hundred iterations while quantum dynamics quantitatively catches it only if h < 3 ■ 10~ 4 
or N > 2 ■ 10 3 whereas 90% derealization is achieved already for h = 0.0009 ( |55|) 
(between figures 15d and 15e). 



5 Summary and conclusions 

In the present paper I have introduced exact unitary quantum Poincare mapping which 
has all the necessary properties: (i) it yields an exact and practically extremely useful 
quantization condition, and (ii) it can be literally interpreted as quantum CSOS-CSOS 
propagator since it is unitary and since energy dependent quantum propagator (Green 
function of the Schrodinger equation) can be decomposed in terms of CSOS-CSOS prop- 
agator and additional three conditional propagators between CS and CSOS which have 
been defined in the paper. In the second part of the paper I have applied this quantum 
SOS method for quantizing a simple but generic 2-dim autonomous Hamiltonian system, 
namely the semiseparable oscillator. I have studied both stationary and time evolving 
(better: SOS evolving) quantum dynamics of the system. Due to extreme efficiency of 
the method, especially for the so-called semiseparable systems, I have been able to go 
orders of magnitude higher than has previously been possible by any other method, up to 
20 millionth eigenstate. Even for geometrically completely generic systems, such as e.g. 
diamagnetic Kepler problem ||I3]| , the method is expected to reach a millionth eigenstate 



25]. I have confirmed Berry- Robnik scenario for level spacing distribution, classification 
of states into regular and irregular class, showed few typical and atypical examples of 
SOS Husimi distributions and wavefunctions of eigenstates, analyzed the distribution of 
zeros of SOS Husimi distributions of eigenstates and found uniformity with cubic near- 
est neighbour repulsion, and successfully compared quantum and classical Poincare SOS 
evolution. 
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Figures 



Figure 1 Localization area distribution V(A) (a), Quantum-classical overlap distribu- 
tion V(B) (b) and contrast distribution V{C) (c) for the 13445 eigenstates of SSO for 
a = 0.03, 6| = 5, bi — 10, ft = 0.0003 (case II). The dark, black curves denote the cumula- 
tive (integrated) distributions, whereas bright, grey curves denote the usual probability 
densities. The thin black curve in (a) denote the cumulative distribution of general- 
ized localization areas A/>(1) which is much less sensitive and therefore less appropriate 
for the classification of eigenstates than the usual localization area distribution (thick 
black curve). The horizontal dotted lines in (a,b) denote the relative volume of classical 
regular component p\ = 0.291 and vertical dotted lines in (a) denote the positions of 
maxima and minima of probability distributions. Scatter diagrams for the corresponding 
three quantities A^,B^,C^, are also shown: vs. (f), A$ vs. C$ (d), and B^ vs. 

(e). Note the two tongues in scatter diagram (f) at small A and at b pa 0.17,0.47 
correspond the regular states (e.g figure 4(g,r),4(h,s)) which live on small islands inside 
chaotic region (figure 4w) and have therefore significant overlap with the chaotic region. 

Figure 2 Same as in figure 1 but for 14231 eigenstates of case I of SSO, a = 0.03, b-\ = 
5, bi — 10, ft, — 0.01. Although P(A) (a) is here still unimodal, P(B) is already bimodal 
and can be used for classification of eigenstates with threshold value of B selected at 
the point where the second (chaotic) peak starts and not at the minimum of probability 
density. 

Figure 3 Cumulative energy level spacing distributions W(S) = Jq dsP(s) for the two 
cases of SSO for a = 0.03, &j = 5, b± = 10. The thick full curve are the numerical 
data, the thin full curve is the best-fit Berry-Robnik distribution, the dashed curve is 
the best-fit Brody distribution, and the dotted curves are the limiting Poisson and GOE 
distributions. For 14231 numerical energy levels of case I (a) ft, — 0.01 one obtains 
significant global fit with Brody distribution, yielding the level repulsion exponent (3 = 
0.142. For 13445 levels in the far semiclassical regime of case II (b) ft = 0.0003 we 
already obtain significant fit by the Berry-Robnik distribution with very accurate value 
of p\ = 0.283. In (c) I show the cumulative level spacing distribution for regular part of 
the spectrum of case II which is indeed very close to Poisson, and (d) for irregular part 
of the spectrum which is also close to GOE (see text for details). 

Figure 4 The figure shows SOS Husimi distributions of six typical eigenstates (above 
the dashed line: a-f,l-q) and five specially picked eigenstates with rare geometric prop- 
erties (below the dashed line: g-k, r-v) for the far semiclassical case II of SSO: a = 
0.03, 6| = 5, b± — 10, ft = 0.0003. In figures (a-k) they are plotted with ten equidistant 
contours while in figures (1-v) they are plotted in logarithmic scale with five contours 
separating regions with different level of greyness, where each two neighbouring contours 
lie by a factor of ?T Q ' 2 = 5.06 apart. The sequential numbers of the eigenstates states 
are around 17684 000 with the following values of the eigenenergies: 0.49999967469 
(a,l), 0.49999969040 (b,m), 0.49999970456 (c,n), 0.49999971683 (d,o), 0.49999974837 
(e,p), 0.49999977984 (f,q), 0.50003435571 (g,r), 0.50000040776 (h,s), 0.49989001264 (i,t), 
0.49996272064 (j,u), 0.49993911720 (k,v). I also give the few (two chaotic and two reg- 
ular) classical SOS orbits in the same scale (w). 
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Figure 5 SOS Husimi distribution of a typical chaotic eigenstate with eigenenergy E = 
0.49999026441 state in a fully chaotic — ergodic regime of SSO (case III: a = 0.25, b[ = 
4, 6| = 11, h — 0.0003). Equidistant (a) and logarithmic (b) contour plots have the same 
parameters as in figure 4. A single overwhelming classical chaotic orbit in SOS of the 
same scale is shown in (c). 

Figure 6 Wavefunctions of a typical chaotic (a) and a regular (b) eigenstate of SSO with 
sequential quantum numbers around 17684 000 (case II: a = 0.03, b^ — 5, — 10, fi = 
0.0003). Their SOS Husimi distributions are shown in figure 4(c,n) for (a) and in figure 
4(d,o) for (b). In the left part of each figure the wavefunction in the entire configuration 
space —1.025 < x < 1.025, — 10 < y < 5 is given with two black windows being magnified 
by a factor of 100 on the upper right side of the figure, in the same left-right order, and 
again with two black windows being magnified by another factor of 100 below, on the 
lower right side of the each figure. The regions where the magnitude of a square of the 
wavefunction (or its average over several tens de Broglie's wavelengths for the left part of 
each figure where the quantum resolution exceeds graphical resolution) is above suitably 
chosen threshold are painted black whereas everything else is white. SOS is indicated 
with a horizontal thin line. 

Figure 7 The wavefunction of a typical chaotic eigenstate with sequential number around 
17684 000 and eigenenergy E = 0.4999026441 (its SOS Husimi distribution is shown in 
figure 5) in a fully chaotic — ergodic regime of SSO (case III: a = 0.25, — 4, b^ — 
11, H = 0.0003). Presentational technique is the same as for figure 6 except for slightly 
extended region of x— coordinate on the left part —1.15 < x < 1.15 due to larger defect. 

Figure 8 The wavefunction of a strongly (bouncing-ball-like) scarred chaotic eigen- 
state with sequential number around 17684 000 and with eigenenergy E = 0.4999999035 
(which is indeed very close (within mean level spacing) to quantization of a lower box 
with 106 nodes in y-direction which gives E = 0.4999999027) of SSO for case II: 
a = 0.03, 6j_ = 5, b-\ = 10, k = 0.0003. Presentational technique is the same as for 
figure 6. 

Figure 9 The zeros of the SOS Husimi distributions are plotted for the same eleven 
eigenstates of SSO in the same scale as in figure 4. The edge of the largest chaotic 
region is marked with a thin curve so that one can easily observe the correlation between 
distribution of zeros and structure of classical SOS. Note the random uniform distribution 
of zeros on the chaotic component and regular distribution of zeros along regular classical 
invariant curves. Note also how zeros avoid regions of maximal Husimi density for regular 
states (a,d,g,h,j) and also for tiny chaotic state (k). 

Figure 10 The zeros of the SOS Husimi distribution for the same typical chaotic eigen- 
state in ergodic regime of SSO as in figure 5 is shown in (a). The superposition of all 
(cca. 34800) zeros of 16 consecutive eigenstates which includes (a) is shown in (b) to 
demonstrate the uniformity of their distribution. For more quantitative conclusions I 
show also (c) the radial density of zeros (in arbitrary units) as determined from the zeros 
of these 16 eigenstates (b). 

Figure 11 The cumulative nearest zero spacing distribution (a) W(S) = J 5 dsZ(s) is 
shown for cca. 16000 "chaotic" zeros of SOS Husimi distributions of 19 chaotic eigen- 
states from the sequence of 24 consecutive eigenstates for SSO of case II: a = 0.03, 6| = 
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5, b i = 10, h = 0.0003, 0.4999995 < E < 0.5 I have also plotted the corresponding T- 
function |26| T(\nS) = ln(— ln(l — W(S))) to analyze behaviour at small and large 
spacings S where T-function is linear for power law and exponential of a power 
exp(— S" 7 ), respectively. The numerical curve (thick curve (a), ±sigma statistical error 
bars (b)) excellently agrees with results obtained from statistical ensemble of GRHD 
(thin curves). The small S and large S parts are magnified in right-lower and left-upper 
window, respectively (a). Dashed line with a slope 4 is a guide for the eye (b) indicating 
that there is a cubic repulsion of zeros at small S. 

Figure 12 The cumulative nearest zero spacing distribution (a) W(S) = dsZ(s) is 
shown for the zeros of SOS Husimi distributions of 16 chaotic states of SSO in the chaotic 
regime (the same eigenstates as in figure 10). Everything else is the same as in figure 11 
including conclusions. 

Figure 13 Quantum Poincare phase space SOS time evolution h n (zf, Zi, E) of an initial 
wave packet located at Xi = 0,p X i = 0.77 is shown (with equidistant contours: a-f, 
and logarithmic contours: g-1 with the same parameters as in figure 4) and compared 
with the corresponding coarse grained classical dynamics f n (zf, E) (equidistant: m-q, 
logarithmic: r-v) for the far semiclassical regime of SSO (case II: a = 0.03, = 5, bi = 
10, h = 0.0003, E = 0.5). 1st (a,g while identical classical initial conditions are not 
plotted), 2nd (b,h,m,r), 3rd (c,i,n,s), 10th (d,j,o,t), 100th (e,k,p,u), and 1000th (f,l,q,v) 
iteration are shown. Note that both, quantum and classical state uniformly condenses 
on the classical invariant chaotic component. 

Figure 14 The same as in figure 13 except for near semiclassical case I of SSO: a = 
0.03, 6 T = 5, fej, = 10, h = 0.01, E = 0.5 and for lst,2nd,3rd,5th,10th, and 100th iteration 
of Quantum Poincare mapping. The quantum state does not condense on the entire 
classical invariant chaotic component of SOS but it remains localized on much smaller 
region of phase space (see figure 15). This phenomenon of quantum localization (in 
autonomous systems) is responsible for fractional power law level repulsion laws (see 
figure 3 a). 

Figure 15 The figure shows localization areas of quantum states which evolve under 
unitary quantum Poincare mapping for the SSO in a diffusively ergodic regime a = 
0.03, 6^ = 500, bi = 1000, E = 0.5 but for eight different values of effective Planck's 
constant: 0.0003 (a), 0.000424 (b), 0.0006 (c), 0.000849 (d), 0.0012 (e), 0.0017 (f), 0.0024 
(g), 0.00339 (h) which increase geometrically by a factor of y/2. The initial state is always 
a wave packet located at Xj = 0,p X i = 0.5. The thick noisy curves denote localization 
areas of quantal states versus the number of iterations while the thick dotted curves 
denote localization areas of coarse grained classical states which evolve under classical 
Poincare mapping. Note the very good agreement between the two curves up to some 
break iteration, where the two curves separate due to quantum localization except in the 
far semiclassical regime (a). The thin (full and dotted) curves denote the same quantities 
but on 10 times smaller iteration scale (up to 33rd iteration). 
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